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ABSTRACT 

It has been suggested that the X-shaped morphology observed in some radio sources can reflect either a recent merger of two su- 
permassive black holes (SMBHs) or the presence of a second active black hole in the galactic nucleus. These scenarios are tested 
by studying the relationship between the black hole mass, radio and optical luminosity, starburst history, and dynamic age of radio 
lobes in a sample of 29 X-shaped radio galaxies drawn from a list of 100 X-shaped radio source candidates identified from the FIRST 
survey. The same relationships are also studied in a control sample consisting of 36 radio-loud active nuclei with similar redshifts 
and optical and radio luminosities. The X-shaped objects are found to have statistically higher black hole masses and older starburst 
activity compared with the objects from the control sample. Implications of these findings are discussed for the black hole merger 
scenario and for the potential presence of active secondary black holes in post-merger galaxies. 

Key words. Galaxies: kinematics and dynamics - Galaxies: formation - Galaxies: nuclei - Black hole physics 



1. Introduction 

X-shaped (or 'winged') radio galaxies are a class of extragalac- 
tic radio sources with two low-surface-brightness radio lobes 
(the 'wings') oriented at an angle to the active, or high-surface- 
brightness, lobes (Leahy & Williams 198 4). Both sets of lobes 
often pass symmetrically through the center of the host galaxy, 
giving the galaxy an X-shaped morphology as seen on radio 
maps. Several authors have attempted to explain the unusual 
structure in X-shaped sources. Rees ( 11978) suggested that the 
X-shaped morphology may result from a precession of the jet 
caused by misalignment between the angular momenta of the 
central supermassive black hole (SMBH) and the gas accreted 
onto it. Subsequently, four different scenarios were proposed for 
the formation of this peculiar radio morphology: 1) backflow 
from the active lobes into the wings (Leahy & Williams 119841 
Capetti et al. |2002); 2) slow conical precession of the jet axis 
(Parma, Ekers & Fanti [19851 Mack et al. [T994t : 3) reorienta- 
tion of the jet axis during which flow continues; and 4) reori- 
entation of the jet axis, but the jet turned off or at greatly re- 
duced power during the change of direction (Dennett-Thorpe et 
al. 120021 1. Merritt & Ekers (120021 1 suggested that the reorien- 
tation of a black hole spin axis owing to a minor merger may 
lead to a sudden flip in the direction of the associated jet. In this 
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case, X-shaped galaxies would be the sites of spin-flips asso- 
ciated with recent coalescences of pairs of supermassive black 
holes (Gopal-Krishna et al. 120031 Gergely & Biermann 120091 . 
The older wings would then represent relics of past radio jets, 
while the active lobes would reflect the activity that ensued after 
the black hole merger. Hence, synchrotron aging should lead to a 
steeper spectral index in the wings or low-surface brightness fea- 
tures than in the high-surface brightness active lobes. However, 
in some X-shaped sources, the wings have a similar or even flat- 
ter spectral index than the active lobes (Lai & Rao l2005"l >. To ex- 
plain this discrepancy Lai & Rao ( 2005) suggested that at least 
some of the X-shaped sources may contain coalescing binary 
SMBH systems with two pairs of jets associated with two unre- 
solved AGN. The activity of the secondary black holes in post- 
merger galaxies is suggested to depend on the separation and 
mass ratio of the two SMBH (Lobanov 2008 ), with equal mass 
mergers providing the longest time during which both SMBH re- 
main active, while minor mergers lead to rapid disruption of the 
accretion disks around the secondary black hole and to quench 
of its activity (the primary black hole in minor merger systems 
may however remain active for a much longer time compared 
with the activity duration in equal mass mergers, cf., Dokuchaev 
[19911) 

Several studies have been carried out in the optical and X- 
ray band, with results in some cases favoring the hydrodynamic 
backflow model over the merger-induced rapid reorientation see- 
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nario (e.g., Hodges-Kluck et al. 120 101 Landt et al. 120101) . and in 
other cases favoring the merger model as the best explanation for 
the observed winged morphology (Hodges-Kluck et al. 12010b ). 
A more detailed review of the different scenarios proposed to 
explain the origin of X-shaped radio galaxies is presented by 
Gopal-Krishna et al. d2QlQb . 

In all schemes in which X-shaped objects are the product 
of galactic mergers, the properties of the nuclear region will be 
affected by the past SMBH merger or the presence of a sec- 
ondary SMBH. This should then be reflected by statistically 
higher masses of the central black holes in galaxies hosting the 
X-shaped radio. If the X-shaped morphology is indeed caused 
by a profound event in the nuclear region, such as a merger of 
two SMBH, then it may also be reflected in the starburst his- 
tory of the host galaxy (Blecha et al. 2010). Investigations of 
the black hole masses and starburst histories can therefore ad- 
dress the question of the physical nature of the X-shaped sources 
and their difference from the "canonical" radio galaxies, if repre- 
sentative samples of both radio types are studied and compared 
against each other. 

In this paper, we determine the black hole masses, lumi- 
nosities, starburst histories, and jet dynamic ages in a sample 
of X-shaped galaxies and compare them with a control sample 
of radio-loud active nuclei with similar redshifts and optical and 
radio luminosities. 

The target and control samples used for the study are intro- 
duced in Sect. 2. Section 3 describes the data analysis and meth- 
ods employed for determination of the black hole masses and 
dynamic ages of the radio emission. Results of these measure- 
ments are presented in Sect. 4. A discussion of the results and 
their broader implication is given in Sect. 5. 

Throughout this paper, we assume a A-CDM cosmology 
with parameters Hq = 73 km s^ 1 Mpc~ l , Q.^ = 0.73 and 
O,,, = 0.27. 

2. The sample 

The sample of AGN analyzed here is drawn from a list of known 
X-shaped sources and a list of 100 'winged' and X-shaped ra- 
dio source candidates (Cheung |20071 i that were morphologically 
identified in the images from the VLA FIRST survey (Becker et 
al. [19951 

Of the 100 X-shaped candidates, 22 were found to have 
spectroscopic information in the Sloan Digital Sky Survey 
(SDSS) data release (DR6; Adelman-McCarthy et al. 120081 . 
and optical spectra for additional 27 objects were obtained by 
Cheung et al. ( 120091 ) with the 9.2 m Hobby-Eberly Telescope 
(HET) at McDonald Observatory and the 6.5 m Multiple-Mirror 
Telescope (MMT) at Mt. Hopkins Observatory. 

Cheung et al. (2009) have defined a subsample of 50 radio 
galaxies with bona fide X-shaped radio morphology. Because 



this selection was based on visual inspection of the candidate 
objects and lacks quantifiable selection criteria, we refer to the 
original sample of Cheung (2007) for selecting appropriate can- 
didates for our study. We select the X-shaped objects with optical 
spectra that show well-detected stellar absorption lines. Our final 
sample comprises 18 X-shaped radio sources with SDSS spec- 
tra complemented by five of the 27 X-shaped objects with spec- 
tra presented in Cheung et al. (120091 ), and six of the previously 
known X-shaped radio sources: 3C192, 4C+32.25, 4C +48.29 
and 1059+169 (Lai & Rao 120071 3C223.1 (Lai & Rao l2005l ) 
and 4C+01.30 (Wang et al.:2003 ) with spectra also available in 
the SDSS. 

In order to evaluate the results obtained for the X-shaped 
radio galaxies, we compile a control sample of 19 radio-loud 
sources from Marchesini et al. ( 120041 ) plus 6 Fanaroff-Riley type 
II (FR II) sources from de Vries et al. (2006), and 11 radio 
loud elliptical galaxies from Gonalez-Serrano & Carballo ( 120001 ) 
that have SDSS spectra, and cover the same ranges of redshift 
(z < 0.3) and optical and radio luminosities as the objects in the 
target sample. The resulting common luminosity ranges for both 
samples are: log AL 5m ^ e [43.0,46.0] and log vL\achz e [39.0, 
44.5]. According to the Kolmogorov-Smirnov test (KS-test) the 
two samples differ at 0.9cr, which warrants making statistical 
comparisons between them. 

3. Data analysis 

The optical spectra of the objects in both samples have been used 
to measure properties of H/3 and [OIII] emission lines and a num- 
ber of stellar absorption lines. These measurements were applied 
to making kinematic mass estimates of the central black holes. 
The absorption line fits were also used to recover starburst his- 
tories in the host galaxies and to estimate the epochs of the most 
recent starburst. Angular sizes of active lobes and inactive wings 
(in X-shaped objects) were measured to determine the dynamic 
ages of the radio emission in the X-shaped objects and in the 
objects from the control sample. 

3. 1 . Stellar absorption lines 

The SDSS spectra collected for the objects in both samples con- 
tain several significant stellar absorption features (such as Ca 
H+K AA3969, 3934, the Mg I b /L15167, 5173, 5184 triplet, 
and the Ca II /L18498, 8542, 8662 triplet, etc.) that can be 
matched against a combination of different synthetic stellar tem- 
plate spectra, yielding an estimate of the stellar velocity disper- 
sion. 

We use the stellar population synthesis code STARLIGHT 
(Asari et al. l2007l Cid Fernandes et al. [2004l[2005l[2007t Mateus 
et al. 120061 ) to model the observed spectrum Ox- The best fit is 
obtained using a linear combination of simple stellar populations 
(SSPs) from the stellar library of Bruzual & Chariot ( 120031 ) and 
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a power-law component representing the AGN continuum emis- 
sion. In our fitting, we apply the standard stellar library consist- 
ing of 150 SSPs, and complement the fit with up to six power- 
law components given by F(A) = 10 20 (A/4020)^, where /3 = -0.5, 
-1, -1.5, -2, -2.5, -3. Each SSP spans across six metallicities, Z 
= 0.005, 0.02, 0.2, 0.4, 1, and 2.5 Z Q , with 25 different ages be- 
tween 1 Myr and 18 Gyr. The Galactic extinction caused by the 
foreground dust screen is modeled and parametrized by the V- 
band extinction, Ay. We adopt the extinction law of Cardelli et 
al. ([19891) . 

The resulting model spectrum M A (combining an SSP and a 
power-law continuum components) is 

N. 

M A (x,Mi ,A v ,v*,o-«) = M Ao C^x j b j , A r A \®G{v t ,cr t ), (1) 

;=i 

where b hA = L s A SP (tj,Zj)/L s A ^ p (tj,Zj) is the f 1 template normal- 
ized at /to, x is the population vector, M Ao is the synthetic flux 
at the normalization wavelength, r A = 10 _0 - 4(A|_Aj o ) is the red- 
dening term, and G(v*,<x„) is the line-of-sight stellar velocity 
distribution, modeled as a Gaussian feature centered at velocity 
v* and broadened by cr». The best fit is reached by minimizing 

x\ 

X 2 (x,M Ao ,A v ,v t ,cr t ) = 2K°-« ~ m a)wa] 2 , (2) 

-1=1 

where the weighted spectrum w A is defined as the inverse of the 
noise in the observed spectra. A more detailed description of the 
STARLIGHT code and its applications can be found in Asari et 
al. d20Q7l >; Cid Fernandes et al. d2004l 120531 [2007b : Mateus et 
al. (120061) ; Leon-Tavares et al. d2QlQb . 

The fit for the X-shaped source J 1424+2637 is shown as an 
example in Fig.Q] The observed spectrum is shown in black, the 
host galaxy model in red, and the AGN power-law continuum in 
blue. The residual spectrum obtained after the subtraction of the 
stellar light and the continuum is shown in green. 

In order to assess the fidelity of the STARLIGHT fit, we in- 
troduce a quality factor Q that combines the reduced^ 2 param- 
eter of the modeled spectra, the velocity dispersion, cr„, and its 
error, 5*: Q = (x 2 8*/cr*y l l 2 . Fits with Q > 10 can be consid- 
ered reliable (these are the fits with reduced^ 2 close to unity and 
fractional errors of the velocity dispersion of < 3 %). 

For some objects with a strong continuum, STARLIGHT can 
fail to synthesize the spectrum of the host galaxy, since the flux 
of the lines of the host galaxy is much fainter than the AGN con- 
tinuum flux. In these cases, we apply the empirical correlation 
cr» = FWHM [OIII]/2.35 obtained by Nelson j2000b assuming 
that cr„ * o- gas . The term FWHM[OIII] describes the full-width- 
at-half-maximum of the narrow component of the [OIII] emis- 
sion line and is determined from a fit to the narrow lines of the 
[OIII]/l/l4959, 5007 A doublet in the residual spectrum. 
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Fig. 1: STARLIGHT fit for the spectrum of the X-shaped source 
J1424+2637. The plot shows the observed spectrum (top black 
line), the modeled spectrum (second red line, displaced from ob- 
served spectrum), the AGN continuum (blue line), and the resid- 
ual spectrum (bottom green line). Note that the residual spec- 
trum contains narrow and (in some objects) broad emission lines. 
Properties of these lines are measured separately (see Sect. 3.2). 



3.2. Emission lines 

We determine the properties of the emission lines using the 
residual spectra obtained after subtracting the host galaxy con- 
tribution and AGN continuum from the STARLIGHT fits. The 
H/? and [OIII] emission lines in these spectra are modeled as a 
combination of the narrow and broad components. We first fit 
three narrow Gaussian components for the oxygen [OIII]4959A 
and [OIII]5007A lines and for the H/? line. Three additional con- 
straints are applied during the fitting to reduce the number of 
free parameters and to provide robust fits: 1) the central wave- 
lengths of the narrow line components are set to their respective 
laboratory wavelengths; 2) the FWHM of the two [OIII] narrow 
components are required to have the same value; 3) the 1:3 ra- 
tio of amplitudes of the [OIII]4959A and [OIII]5007A narrow 
components is fixed. 

After adjusting the narrow components, broad Gaussian 
components are added to the H/3 line, and the joint fit is further 
adjusted until the relative residuals are reduced to below 0.02- 
0.03. 

3.3. Starburst histories 

The STARLIGHT model for the observed spectra also yields the 
light fraction xj, mass fraction Minij, age tj, and metallicity Zj, 
of the stellar populations used in the fit. We use these parame- 
ters to derive starburst histories and estimate the epochs of the 
most recent starbursts in the studied galaxies. We apply Gaussian 
smoothing to the individual starburst events (see Fig.|2]i and de- 
termine the epoch of the most recent starburst episode. 
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Fig. 2: STARLIGHT fits for J1327-0203: the light fraction (jcj) 
vs. age (top) and mass fraction Minij vs. age (bottom) for each 
stellar population used. The squares correspond to stellar pop- 
ulations with different metallicities (in the color version: blue 
0.005 Z Q , purple 0.020 Z Q , green 0.200 Z G , yellow 0.400 Z Q , 
orange 1.000 Z and red 2.500 Z . The solid curve represents 
Gaussian smoothing applied to the mass and light fraction distri- 
butions of individual starbursts. The stars indicate the identified 
resulting starburst episodes. 

3.4. Optical continuum 

Most of the X-shaped radio sources and some objects from the 
control sample are type II AGN (weakly beamed radio sources), 
for which STARLIGHT cannot provide reliable estimates of the 
AGN continuum flux. We estimate the rest-frame continuum flux 
at 5100 A from the SDSS photometry, with the relation (Wu & 
Liu l2004T > 

-te-r)/2.5/o g (623 1/4770) 

/WJy] = 3631xl0- u ^ rAnn ~~. ■ (3) 



5100(1 + z) 

where z is the redshift and the g and r fiber magnitudes are ob- 
tained from the SDSS and corrected for the Galactic extinction 
Ay (taken from Schlegel et al. 11998b . It should be noted that the 
flux obtained with this method also contains a contribution from 
the host galaxy, but this does not introduce a strong bias in our 
estimates. 

To assess the spectral classification of the X-shaped and con- 
trol host galaxies, we measure the Ca II break of their absop- 
tion optical spectrum. This break is typically seen in the spec- 
trum of elliptical galaxies and is described by a factor Ccaii = 
(/V - /-)//+, where f- and /+ are the fluxes in the rest frame 
wavelength regions 3750-3950 A and 4050-4250 A, respectively 



(Landt et al. l2002l >. The Ca II break was used as an additional 
criterion to separate blazars from radio galaxies. Stocke et al. 
(119911) adopted a maximum value of Ccaii = 0.25 for BL Lacs 
to ensure the presence of a substantial non-thermal jet continuum 
in addition to the thermal spectrum of the elliptical host galaxy. 
This limit was later increased to Ccaii < 0.4 by other authors 
(Marcha et al. [19961 Plotkin et al. 120051 . 

Landt et al. (120101) provide the Ca II break values of 16 of 
the X-shaped sources studied here. For the rest of the X-shaped 
sources and control objects, we measure the C factor of the rest- 
frame optical spectra using the IRAF task guiapps.spectool. 

3.5. Black hole masses 

The measured stellar velocity dispersion, cr„, can be connected 
with the mass, Mbh, of the central black hole through an empir- 
ical relation (Gebhardt et al. 120001 Tremaine et al. 120021 ): 



M BH = 1.349 xl0 8 M c 



f \ 

\200kms~ 1 ; 



4.02±0.32 



(4) 



. This relation is valid under the assumption that the kinematics 
of the stars in the bulge of the host galaxy is dominated by the 
gravitational potential of the central SMBH (Ferrarese & Merritt 
1200Tb . 

3.6. Dynamic age of radio lobes 

The dynamic age, f a , of the high-surface-brightness (active) ra- 
dio lobes is obtained from their angular size, 6> a , which is de- 
fined as the separation between the center of the radio source 
and the most distant contour in the FIRST image. The dynamic 
age is then given by 6 a /v a , where v a is the lobe advance speed 
(we adopt the commonly assumed v a « 0.1 c; cf., Tingay et al. 
119981 ). Assuming that the fuelling of the low-surface-brightness 
lobes of the X-shaped sources had stopped after the high-surface 
ones were activated, the dynamic age of the passive lobes t p dur- 
ing their active stage can be estimated as 



(5) 



where 6 p is the angular size of the low-surface brightness lobes 
and v p is their expansion speed during the inactive stage. In the 
absence of an observational estimate for v p , we use v p = 0.01 c in 
our calculation. It should be noted that reducing v p further only 
has a small effect on the ages derived for the passive lobes. The 
total dynamic age of the X-shaped sources can then be obtained 
from the sum of the dynamic ages of the active and passive lobes. 

4. Results 

The combined results from fitting the optical spectra, black hole 
mass calculations, and age measurements for the radio lobes and 
most recent starbursts are presented in Tables 2-3. 
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Tables 2 and 3 list (for the X-shaped objects and the control 
sample, respectively) the object name based on J2000.0 coor- 
dinates (Col. 1), other common catalog names (Col. 2), stellar 
velocity dispersion (Col. 3), black hole mass derived from <x, 
(Col. 4), optical luminosity of the AGN (Col. 5) and of the host 
galaxy (Col. 6), radio luminosity (Col. 7), dynamic age of the 
radio lobes (Col. 8), most recent starburst age (Col. 9), spec- 
troscopic redshift (Col. 10), quality factor of the STARLIGHT 
fit (Col. 1 1), the value of Ca II break factor (Col. 12), and refer- 
ences for the control sources (Table 3, Col. 13). For the X-shaped 
sources, the total (active + passive lobe) age of the radio emis- 
sion is given in brackets in Col. 8. 

The control sample objects 1217+023,2349-014, 1004+130, 
3C277.1 and 3C254 are quasars with strong power-law contin- 
uum in their spectra. Estimates of cr„ for these objects are ob- 
tained from the measured width of the [OIII] line as described 
in Sect. 3.1. The STARLIGHT results for two X-shaped sources 
(J0941-0143 and J1348+4411) have provided unreliable cr, fits 
(with values values of cr* below the instrumental resolution of 
70km/s in the SDSS spectra, and quality factors Q <10) and are 
excluded from further analysis. 
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Fig. 3: Optical continuum luminosity versus radio luminosity at 
1.4 GHz for X-shaped (crosses) and control sources (squares). 
A small sub-region ("Region 1" marked by a rectangle inside 
the plot) is identified to provide a tighter luminosity match be- 
tween the two samples. In the further analysis, Region 1 is also 
subdivided into two equal bins in the radio luminosity to assess 
a possible luminosity dependence of the black hole mass esti- 
mates. 



4.1. Luminosity matching 

Luminosity matching between the target and control samples is 
illustrated in Fig. 3, where the continuum luminosities derived 
from the SDSS magnitudes (Eq.[3) are plotted against the radio 
luminosities at 1 .4 GHz. The entire original radio-optical lumi- 
nosity range, including all sources, is called Region hereafter. 
To provide a tighter match between the samples (at the expense 
of reducing the sample sizes), we define a smaller subregion or 
window (shown in Fig. [3} called Region 1. Region 1 ranges from 
log^L 510oA e [43.5, 44.25] to log vL lMHz e [40.25, 42.5] and 
mainly excludes the control sources with the lowest radio lu- 
minosities and highest optical luminosities because there are no 
X-shaped sources with such luminosities. The KS-test indicates 
that the statistical difference between the two samples in Region 
1 is 0.9<x for the optical luminosity and 1 .Oct for the radio lumi- 
nosity. 

4.2. Spectral classification 

We base our spectral classification on color-color separation and 
on the Ca II break measurements in the target and control sam- 
ples. All X-shaped sources analyzed have Cc a n > 0.25, with 
92% of the sources (23 out of 25) having Cc a n — 0.4. These 
values ensure that the host galaxy of the X-shaped objects is el- 
liptical and that it is dominated by the thermal spectrum, without 
any significant contribution of the non-thermal jet continuum. 

For the control sample, however, the Ca II break can be de- 
termined for only 20 out of the 36 sources. Ninety percent of 
these sources (18 out of 20) have C Ca n > 0.25, and 70% (14 out 



of 20) have Cc a n — 0.4. Only 56% of the control sources there- 
fore have a thermal spectrum that dominates the host elliptical 
galaxy. 

The spectral classification of the X-shaped and control 
sources can be further refined with color-color diagrams. We use 
the g, r and u magnitudes obtained from SDSS to plot a g-r vs u- 
g color-color diagram (Fig.|4]i for the two samples. We include in 
this plot only the sources with a spectrum dominated by the host 
galaxy (as determined from the Ca II break). According to the 
distribution of galaxies in the u-g vs g-r diagram from Strateva 
et al. (120011 1. all galaxies lying above the u - r — 2.22 separator 
line are elliptical systems. This is the case for all the X-shaped 
radio sources and for all the control sources for which the SDSS 
magnitudes are available, with only one source from each of the 
samples lying slightly below the u — r = 2.22 separator line. 
The X-shaped source J0049+0059 (with u - g = 5.0) is not plot- 
ted in Fig.|U but it has been classified as an elliptical galaxy as 
well. This spectral classification agrees with the results obtained 
using the Ca II break feature and accordingly confirms that all 
X-shaped analyzed here are hosted by elliptical galaxies. 

In order to better study the relation between the galactic hosts 
of the X-shaped and control radio sources, we will also consider 
a control subsample of elliptical galaxies from now on. This sub- 
sample is defined according to the galactic type obtained by the 
SDSS color separation described above and the Ca II break, and 
contains the 20 control sources for which the Ccaii value could 
be determined. 
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Fig. 4: Color-color diagram (g-r colors versus u-g colors) for X- 
shaped sources (cross) and control sample (square). Black line: 
u — r = 2.22 galaxy type separator from Strateva et al. (120011) . 
Sources situated above the line are classified as ellipticals. 



4.3. Black hole masses 

In Table 1 we show the mean black hole mass for each sam- 
ple and subsamples and for each of the regions described above 
(Regions and 1). The ratios of the mean black hole mass de- 
fined as r xc = <M B H,X-shaped>/<^BH .control 

> are given as well. 

A comparison of the black hole masses derived for the entire 
samples (objects in Region 0) yields ratios of mean BH mass 
ranging between 1.48 for X-shaped/control (all) to 1.13 for X- 
shaped/control (ellipticals). 

The mass ratios derived for objects in Region 1 (tighter com- 
mon range of luminosities) are higher, with the mean BH mass 
of the X-shaped sample being nearly twice the one derived for 
the whole control sample. In order to further reduce any potential 
effect of the source luminosity on the derived black hole mass, 
we calculate mass ratios in two equal-size bins in Region 1 : bin 
1 ranges from log vLi.4c//: 6 [40.25, 41.4] and bin 2 from log 
vLiaghz e [41.5, 42.5]. One can see from the resulting ratios 
given in Table 1 that the mass ratio does not show any evolution 
between the bins for the whole control sample. The decrease of 
the mass ratio in bin 2 for the subsample of control ellipticals is 
owing to the small number of sources in this bin (3 sources). 

The differences in the black hole masses are further illus- 
trated in Figs.|3H6l which show histograms of the mass distribu- 
tions for X-shaped and control samples in Regions and 1, re- 
spectively. In both regions the X-shaped objects show a tendency 
to have higher black hole masses. This trend is more significant 
in the tighter Region 1, where 60% of the X-shaped sources have 
log Mbh > 8.25 M while only 25% of the whole control sample 
does. . 

The KS-test applied to the BH masses of the X-shaped and 
control objects indicates that the mass distributions are different 
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Fig. 5: Histogram of the black hole mass in Region for X- 
shaped sources (top), control all (middle), and control ellipticals 
(bottom). 
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Fig. 6: Histogram of the black hole mass in Region 1 for X- 
shaped sources (top), control all (middle), and control ellipticals 
(bottom). 



at a statistical significance of 1 .6<x. The differences are slightly 
larger for objects in Region 1 (1.9cr). 

The statistically larger black hole mass of the X-shaped sam- 
ple implies that these objects are possibly located in galaxies that 
have undergone strong major activity in the past, with either one 
major merger event or multiple minor mergers. 

4.4. Starbursts 

The starburst and radio lobe ages can be used to see whether 
starburst activity is different in the X-shaped objects and whether 
it can be related to the production of radio emission in both the 
target and control samples. 

A histogram of the ages of the most recent starburst (Fig. [7]) 
in the objects from Region shows that a sizeable fraction of ob- 
jects in both samples have relatively recent starbursts, with ages 



M. Mezcua et al.: Black hole masses and starbursts in X-shaped radio sources 
Table 1 : Mean black hole masses and mass ratios 



7 



Region 



X-shaped 



Control (all) Control (ellipticals) 



(M BH )[M ) 



21.81+^2 x io 7 



14.78+^5 x 10 7 



1 48 +0 - 24 



19.38^ x 10 7 



i i 1+0.20 
i - iJ -().17 



(M BH )[M ] 



21.03^ x 10 7 1 0.87+J-34 x 10 7 



i Q4+0.42 



13.98+^8 x 10 7 



1 50 +(I36 
i .j»_ 29 



1 Bin 1 



Bin 2 



(M BH )[M ) 



(M BH )[M G ) 



23.39^-g x 10 7 



15.05+^3 x 10 7 



i cc+0.44 
l.JJ_ 034 



16.00^-2" x 10 7 8.75+|^ x 10 7 

1 ot+0.57 

1 - 8 - ) -n.43 



14.10^ °2 x 10 7 



1 66 +a59 



14.03+|;^ x 10 7 



I 14+0.72 



of 10 60 -10 years, which are possibly related to the current 
jet activity. The distribution of starburst ages in the X-shaped 
sample is much broader, with half of the objects having star- 
bursts older than 10 8 years, while only 19% of all the control 
sources and 25% of the control ellipticals exhibit these star- 
bursts. Because the largest dynamic age of the active lobes is 
~ 10 7 years (see Tables 2 and 3, Col. 8), any starburst activity 
significantly older than 10 7 years is not likely to be related to 
the active lobes. It can be speculated that these older starbursts 
are related to galactic mergers themselves or to the putative coa- 
lescence of the central black holes in post-merger galaxies. The 
latter possibility can be viably tested by comparing the starburst 
ages to the dynamic ages of the active and passive lobes in the 
X-shaped objects. 

The KS-test applied to the distributions of the most recent 
starburst ages gives a statistical significance of 2.3cr for the two 
samples being different. To account for a possible dependence of 
this difference on the galactic type, we apply the KS-test to the 
X-shaped sample and the subsample of control ellipticals. The 
results of the test show that the starburst ages of the X-shaped 
sources and the control ellipticals are still different at a statistical 
significance of 2.1cr. 

In Fig. [8] histograms of the logarithm of the ratio of the 
dynamic age and most recent starburst age are plotted for the 
X-shaped sources and control sources of Region 0. The mean 
logarithmic ratios are -1.29 + 0.23 and -0.14 + 0.18 for the 
X-shaped objects and the control sample, respectively. The X- 
shaped sources tend to have starburst ages that are older than 
the dynamic ages of the radio lobes, while in both the control 
sample and the control subsample of ellipticals these ages are 
comparable. The starburst activity in X-shaped sources is there- 
fore likely not related to the active lobes. The KS-test shows 
that the ratio distributions are different at a statistical signifi- 
cance of 2.8<x. This difference in the starburst/dynamic age ra- 



tios may support the scenario in which the active lobes of the 
X-shaped sources are due to a possible reorientation caused by 
a black hole merger (Merritt & Ekers 2002) that leaves the old 
low-surface-brightness lobes inactive. Assuming that the low- 
surface-brightness lobes became inactive when the high-surface 
ones were activated, the dynamic age of the passive lobes during 
their active stage can be determined using Eq.[5] The ratio of the 
total dynamic age of the active plus passive lobes to the starburst 
age is plotted in Fig. [8] (second panel), and it indicates that the 
starburst age still remains older than the total dynamic age of 
the lobes. This suggests that the starburst activity in X-shaped 
sources had occured before the possible reorientation owing to a 
black hole merger, and it may have been related to the galactic 
merger itself. 

Instead of analyzing all the X-shaped sources with optical 
spectra available, we can constrain the X-shaped sample to only 
the X-shaped radio galaxies included in the bona fide sample 
of Landt et al. (|20T0T >. This implies taking out 10 of the 29 X- 
shaped sources included in our sample, which leads to a 50% 
increase in the statistical erros, but does not change the obtained 
results. The statistical studies would be improved, on the other 
hand, by the addition of more X-shaped radio sources to the sam- 
ple. This will be done with the availability of new optical spectra. 

5. Summary 

We determined and studied the black hole masses and starburst 
histories of 29 X-shaped radio sources to attempt to find out the 
most likely physical scenario for the formation of their peculiar 
morphology. We compared the black hole masses and proper- 
ties of the starburst activity derived for this sample with those 
obtained for a control sample of 36 galaxies with common opti- 
cal and radio luminosities. The main results of this study can be 
summarized as follows. 
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Fig. 7: Histogram of the most recent starburst for X-shaped 
sources (top), control all (middle), and control ellipticals (bot- 
tom). 
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Fig. 8: Logarithmic ratios of the dynamic age of the radio emis- 
sion to the age of the most recent starburst for Region 0. The top 
row presents the ages of active lobes in X-shaped objects, while 
the second row shows the sum of the ages of the active and pas- 
sive lobes. The two bottom rows show distributions of the age 
ratios in the control sample (third row) and the subsample of 
control ellipticals (fourth row). 



The Ca II break feature of the absorption spectrum and 
the optical colors obtained from the SDSS u,g, r photome- 
try data indicate that all X-shaped radio sources studied in 
this paper are hosted by elliptical galaxies. In the hierarchi- 
cal galaxy-formation model, today's galaxies are the product of 
frequent galaxy merging, triggering starburst and active galac- 
tic nuclei activity and possibly forming supermassive binary 
black hole systems (Frenk et al. 119881 Hutchings & Campbell 
[19831 Sanders et al. [19881 Hopkins et al. l2006l l. There appears 
to be supporting evidence that elliptical galaxies have formed 
from mergers in a hierarchical way (Toomre & Toomre 119721 
Efstathiou & Silk 1983). Therefore the trend for X-shaped ra- 



dio sources to reside in elliptical galaxies suggests that their X- 
shaped morphology could be related to the galaxy merger sce- 
nario. 

The mean black hole mass of the objects in the X-shaped 
sample is found to be statistically higher than in the control sam- 
ple. This result holds for the entire range of the radio and optical 
luminosities as well as for smaller subranges and for exact galac- 
tic type-matching between the two samples. 

The statistically higher black hole masses obtained for the 
X-shaped objects suggest that their peculiar radio morphology 
can indeed result from a reorientation of the jet axis owing to 
black hole coalescence or from the presence of two active central 
engines in the nuclear region. 

A comparison of the starburst ages and dynamic ages of the 
radio lobes reveals that the most recent starburst activity in X- 
shaped sources is older than in the control sample and that X- 
shaped radio sources had their most recent starburst before their 
active lobes were formed. These results lend further support to 
the scenario in which the high-surface-brightness lobes might 
have become active due to reorientation caused by a black hole 
coalescence, while the peak of starburst activity could be con- 
nected to the galactic merger itself. 

Another scenario proposed to explain the X-shaped mor- 
phology in X-shaped radio galaxies employs the backflow emis- 
sion of the active lobes into the wings. This model is based on the 
presence of pressure gradients in the confining medium of the jet 
and does not seem to be able to explain the results we obtained. 
Although our findings argue in favor of the merger scenario, they 
do not rule out entirely the backflow model. It should be further 
noted that it is quite possible that some X-shaped sources may 
indeed have undergone a merger, while in others the peculiar 
shape might be caused by a backflow of the active jets. 
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Table 2: X-shaped objects 



Name 




log M BH 


log 


log u£ 


log vL rad 


log t. 


. (log fp) 


log t s b 


; 


Q 


Ccall 


T7000 




[knVs] 


[M e ] 


[erg/s] 


[erg/s] 


[erg/s] 




[yr] 


[yr] 










(7\ 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


(12) 


toooi oon 




277.32 ± 4.55 


8.70 ± 0.03 


-b 


43.77 


41.27 


6.37 


(6.60) 


6.50 


0.247 


83.94 


0.43' 


Tnn/iQ_i_nn^o 




282.23 ± 9.64 


8.73 ± 0.06 


-b 


43.43 


41.84 


6.68 


(6.83) 


9.20 


0.304 


31.07 




Tncnc_i_0/ino 
JUoUo+Z4UV 




209.94 ± 7.97 


8.21 ± 0.07 


-b 


43.85 


41.82 


6.52 


(6.71) 


6.70 


0.060 


41.46 


0.36 


TfiS 1 ^_i_A^A7 




251.25 ± 5.79 


8.53 ± 0.04 


-b 


44.11 


41.31 


6.29 


(6.49) 


9.10 


0.128 


44.69 


0.45^ 


TflG3 1 1 O 




194.80 ± 6.21 


8.08 ± 0.06 


42.42 ± 0.14 


43.97 


41.26 


6.69 


(6.85) 


6.50 


0.051 


59.44 


0.44 / 


TAQIO . QOCO 

JUojo+JZjJ 




241.60 ± 11.43 


8.46 ± 0.08 


-b 


43.92 


41.37 


6.59 


(6.76) 


8.50 


0.213 


28.77 


0.42 






250.37 ± 12.90 


8.52 ± 0.09 


42.33 ± 2.65 


— C 


42.13 


6.08 


(6.65) 


8.00 


0.356 


24.55 




TOQ7zl-i-,d7H 




231.16 ± 15.72 


8.38 ±0.12 


-b _ 


43.58 


41.84 


6.54 


(6.77) 


9.00 


0.227 


15.96 


0.42' 


toqai oiai 




154.92 ± 20.61 


7.68 ± 0.23 


42.48 ± 5.15 


43.19 


42.84 


6.14 


(6.64) 


7.00 


0.384 


4.92 e 




TOQA1 o-^QAJ. 
J Uyt Itj ^H-H- 


^P993 1 
jLZZj.l 


196.90 + 9.14 


8.10 ±0.08 


42.30 ± 0.87 


43.94 


41.88 


6.44 


(6.78) 


9.11 


0.107 


35.30 


0.38 / 


ti oo^-u 1 1 




272.88 ± 6.46 


8.67 ± 0.04 


-b _ 


44.12 


41.37 


6.43 


(6.67) 


9.10 


0.166 


66.08 


0.48^ 


t i ooo-lAsi 1 




194.30 ± 7.94 


8.08 ± 0.07 


42.28 ± 0.19 


43.78 


41.22 


6.72 


(7.02) 


8.46 


0.053 


38.82 


0.44 / 


t i oao-l^o^ 

J 1U4Ut jUjO 




211.57 ± 15.88 


8.23 ±0.13 


-b _ 


43.77 


41.24 


6.44 


(6.65) 


9.20 


0.154 


20.74 


0.42 


T 10414-1111 




195.27 ± 8.24 


8.09 ± 0.07 


41.15 ± 1.15 


43.63 


40.49 


5.80 


(6.09) 


6.80 


0.036 


32.66 


0.41 


T 1 101-1-1 f\df) 


1UJ7T 1U7 


220.98 ± 6.42 


8.30 ± 0.05 


41.01 ± 5.18 


43.69 


41.02 


6.62 


(6.98) 


6.00 


0.069 


60.03 


0.45^ 


Ti 1 1 1 




254.75 ± 4.96 


8.55 ± 0.03 


-b _ 


44.17 


41.19 


6.16 


(6.43) 


6.50 


0.074 


74.31 


0.43 


TI 1 10-uOOSR 
J 1 IjUtUUjo 


An _i.ni 


151.86 + 6.78 


7.65 ± 0.08 


43.12 ± 0.15 


44.10 


41.65 


6.24 


(6.61) 


9.40 


0.132 


24.88 


0.26 / 


T1 1 AO-l 1 0^7 
J 1 14Ut IUj / 




196.37 ± 7.80 


8.10 ±0.07 


40.51 ±20.11 


43.61 


41.09 


6.28 


(6.51) 


6.00 


0.081 


40.99 


0.4^ 


T 1 707-i_11^7 
J 1ZU / tjjjZ 




181.28 ± 3.79 


7.96 ± 0.04 


41.85 ±0.84 


44.13 


41.02 


6.19 


(6.44) 


6.00 


0.079 


51.36 


0.28^ 


J 1Z1U— UJ41 




221.29 ± 9.90 


8.30 ± 0.08 


42.14 ±2.94 


43.90 


41.62 


5.97 


(6.56) 


8.70 


0.178 


33.61 


0.40^ 


J 1Z1U+1 1Z1 




223.51 ±9.41 


8.32 ± 0.07 


-b _ 


43.96 


41.35 


6.66 


(6.78) 


7.60 


0.196 


33.55 


0.42 


ti i?7-o?oi 




237.30 ± 9.27 


8.43 ± 0.07 


-b _ 


44.01 


42.22 


6.39 


(6.65) 


7.70 


0.183 


38.98 


0.48 7 


J1330-0206 




219.35 ± 7.72 


8.29 ± 0.06 


-b _ 


44.04 


40.59 


6.35 


(6.61) 


7.00 


0.087 


42.40 


0.40 


J1339-0016 




323.77 + 4.55 


8.97 ± 0.02 


-b _ 


44.15 


41.22 


6.55 


(6.81) 


8.50 


0.145 


88.09 


0.44 


J1348+4411 




108.42 ± 27.71 


7.05 ± 0.45 


41.67 ± 10.43 


43.28 


41.63 


5.87 


(6.53) 


6.90 


0.267 


2.95 e 




J 1424+2637 




174.99 ± 6.53 


7.90 ± 0.07 


40.48 ± 6.37 


43.97 


40.58 


6.13 


(6.37) 


6.50 


0.037 


45.36 


0.40 


J1444+4147 




226.06 ± 9.35 


8.34 ± 0.07 


-b _ 


43.74 


41.74 


6.60 


(6.80) 


8.50 


0.188 


35.89 


0.44^ 


J1455+3237 




224.17 ± 5.52 


8.33 ± 0.04 


41.57 ± 1.78 


43.84 


40.42 


6.10 


(6.31) 


6.50 


0.084 


68.39 


0.39 


J1614+2817 




344.13 ± 6.21 


9.08 ± 0.03 


41.32 ±2.65 


44.02 


41.24 


5.59 


(6.22) 


6.00 


0.108 


25.85 


0.42 / 



Column designation: (1) - object name based on J2000.0 coordinates; (2) - other common catalog names; (3) - stellar velocity dispersion obtained from STARLIGHT; (4) - black hole mass 
obtained from cr,; (5) - 5100 A continuum luminosity from STARLIGHT; (6) - 5100 A continuum luminosity from SDSS photometry; (7) - 1.4 GHz radio luminosity; (8) - dynamic age of the 
active (active+passive) lobes; (9) - age of the most recent starburst; (10) - spectroscopic redshift from SDSS; (11) - quality factor; (12) value of Ca II break factor . Notes: a - velocity dispersion 
obtained using the correlation cr, =FWHM [OIII]/2.35.; b - STARLIGHT could not fit the continuum luminosity; c - no SDSS photometry available; d - no STARLIGHT fit; e - low fidelity of the 
STARLIGHT fit; / - values of Ca II break from Landt et al. ( 120101 
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GOO 


J092 1+4538 


3C219 


188.98 ± 8.49 


8.03 ± 0.08 


42.52 ± 0.07 


43.96 


43.07 


6.94 


6.00 


0.174 


21.14 


- 


M04 


J0925+3627 


4C36.14 


256.37 ± 0.35 


8.56 ± 0.04 


— b 


44.26 


39.36 


6.60 


7.00 


0.112 


68.51 


0.43 


GOO 


J0939+3553 


3C223 


181.42 ± 8.26 


7.96 ± 0.08 


42.45 ± 0.06 


43.81 


42.47 


7.06 


9.40 


0.137 


24.50 


- 


M04 


J094 1+5751 


J094 1+5751 


189.32 ± 7.83 


8.03 ± 0.07 


— b 


43.97 


40.46 


6.43 


6.00 


0.159 


29.56 


- 


V06 


J0947+0725 


3C227 


112.51 ± 3.35 


7.13 ± 0.05 


43.63 ± 0.00 


43.95 


41.47 


6.78 


6.00 


0.086 


24.43 


- 


M04 


J 1006+2554 


B2 1003+26 


240.82 ± 0.40 


8.45 ± 0.06 




44.52 


40.52 


7.17 


9.11 


0.117 


38.18 


0.44 


GOO 


J 1006+3454 


3C236 


247.41 ± 7.26 


8.5 ± 0.05 




44.17 


42.07 


7.74 


6.50 


0.101 


54.96 


0.33 


M04 


J1007+1248 


1004+130 


167.82° ± 1.05 


7.77 ± 0.01 




45.46 


42.49 


6.89 




0.241 




- 


M04 


J1031+5225 


J1031+5225 


204.85 ± 11.06 


8.17 ± 0.09 


42.34 ± 0.10 


43.69 


42.01 


6.42 


6.00 


0.166 


22.09 


0.28 


V06 


J 1040+2957 


4C30.19 


161.97 ± 0.06 


7.76 ± 0.04 


—b 


44.11 


41.03 




8.01 


0.091 


61.46 


0.18 


GOO 


J1055+5202 


J 1055 +5202 


162.90 ± 12.40 


7.77 ± 0.13 


42.61 ±0.11 


43.78 


41.74 


6.46 


6.00 


0.187 


12.12 




V06 


Jl 105+3009 


B2 1102+3OA 


ZOZ.ZO + U.33 


5.0U ± U.U4 


A 1 1 A _i_ 


A A 1^ 
44.33 


1Q 11 

3y.zz 


0.33 


0.3U 


U.U/2 


HA 11 

/4.2Z 


n A 1 
U.43 


uuu 


Jl 114+4037 


3C254 


11A IfV! _i_ 1 TO 


8.35 ± 0.03 


-b 


45.81 


44.3 1 


1.55 


-d 


0.736 


-d 




M04 


Jl 154+0238 


Jl 154+0238 


169.59 ± 11.91 


/.84 ± 0.12 


-b 


43.66 


40.64 


6.46 


6.70 


0.21 1 


18.78 


0.31 


V06 


J1219+0549 


3C270 


269.46 ± 3.72 


8.65 ± 0.02 


-h 


44.04 


40.63 


6.12 


7.00 


0.007 


84.46 


0.41 


M04 


J1220+0203 


1217+023 


174.94" ± 1.72 


7,84 ± 0.02 


-b 


45.25 


42.16 


6.84 


-d 


0.240 


-d 




M04 


J1252+5634 


3C277.1 


235.90" ± 0.57 


8.37 ± 0.00 


-b 


44.62 


43.10 


5.04 


-d 


0.320 


-d 


- 


M04 


J1259+2757 


NGC4874 


251.47 ± 0.32 


8.53 ± 0.04 


—h 


44.36 


39.59 


6.30 


6.00 


0.024 


57.54 


0.42 


GOO 


J1319+2938 


4C29.47 


198.73 ± 0.17 


8.12 ± 0.06 


-b 


44.09 


41.45 


6.35 


8.01 


0.073 


53.97 


0.39 


GOO 


J 132 1+4235 


3C285 


162.58 ± 9.46 


7.77 ±0.10 


—b 


43.86 


41.68 


6.6 


6.00 


0.079 


26.37 


0.24 


M04 


J 1332+0200 


3C287.1 


246.02 ± 7.86 


8.49 ± 0.06 


43.24 ± 0.04 


44.01 


40.27 


6.82 


6.00 


0.216 


32.12 


- 


M04 


J1341+5344 


J1341+5344 


225.84 ± 7.07 


8.34 ± 0.05 


—b 


43.87 


40.02 


6.55 


6.00 


0.141 


42.80 


0.31 


V06 


J1350+2816 


B2 1347+28 


209.58 ± 0.23 


8.21 ± 0.06 


-b 


43.98 


40.65 


6.34 


8.01 


0.072 


46.69 


0.41 


GOO 


J1430+5201 


3C303 


167.41 ± 9.77 


7.82 ±0.10 


43.22 ± 0.01 


43.96 


42.32 


6.28 


7.20 


0.141 


16.86 


- 


M04 


J1512+0203 


J1512+0203 


194.50 ± 10.45 


8.08 ± 0.09 


42.79 ± 0.09 


44.01 


42.30 


6.52 


6.00 


0.219 


25.33 




V06 


J1529+3042 


1527+30 


322.89 ± 7.49 


8.97 ± 0.01 


-b 


44.29 


40.12 


6.38 


8.20 


0.114 


64.60 


0.46 


GOO 


J1559+2556 


B2 1557+26 


223.51 ±0.24 


8.32 ± 0.05 


41.09 ± 15.68 


43.93 


39.89 




6.00 


0.045 


54.67 


0.4 


GOO 


J1611+3103 


B2 1609+31 


199.34 ± 0.23 


8.12 ±0.08 


-b 


43.93 


40.70 


5.74 


6.00 


0.095 


31.15 


0.42 


GOO 


J1615+2726 


1613+27 


221.13 + 6.13 


8.31 ±0.05 


41.11 ± 0.40 


43.99 


40.35 


5.89 


6.00 


0.065 


58.17 


0.42 


M04 


J1617+3222 


3C332 


172.37 ± 3.80 


8.64 ± 0.01 


40.90 ± 4.45 


44.05 


42.41 


6.62 


6.00 


0.151 


33.13 




M04 


J2351-0109 


2349-014 


201.97" ±2.12 


8.72 ± 0.01 


-b 


44.95 


42.28 


6.06 


-d 


0.174 


-d 




M04 



Column designation: (1) - object name based on J2000.0 coordinates; (2) - other common catalog names; (3) ; - stellar velocity dispersion obtained from STARLIGHT; (4) - black hole mass obtained from cr,; (5) - 5100 
A continuum luminosity from STARLIGHT; (6) - 5 100 A continuum luminosity from SDSS photometry; (7) - 1 .4 GHz radio luminosity; (8) - dynamic age of the active lobes; (9) - age of the most recent starburst; (10) - 
spectroscopic redshift from SDSS; (1 1) - quality factor; (12) value of Ca II break factor; (13) - Original references. Notes: a - velocity dispersion obtained using the correlation o",=FWHM [OIII]/2.35.; b - STARLIGHT 
could not fit the continuum luminosity; c - not available; d - no STARLIGHT fit. REFERENCES.- (GOO) Gonalez-Serrano & Carballo <2"0TjQl ; (M04) Marchesini et al. <2"0Tj4l ; (V06) de Vries et al. (20061 



